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>. 

^O . We propose a simple extension of the MSSM based on extra compact dimensions 

which includes an SU{2) (^ U{1) singlet superfield. The fermion present in this 

Q^^ ' superfield is the sterile neutrino, which combines with one linear combination of 

^^ ■ f e — f ^ — z^r to form a Dirac pair whose mass accounts for the LSND anomaly. Its 

small mass can be ascribed to a volume suppression factor associated with extra 

compact dimensions. On the other hand the sterile neutrino scalar partner can 

""q,' trigger the spontaneous violation of R-parity, thereby inducing the necessary mass 

' , splittings to fit also the solar and atmospheric neutrino data. Thus the model can 

g_) \ explain all neutrino oscillation data. It leads to four predictions for the neutrino 

>-Ch ■ oscillation parameters and implies that the atmospheric neutrino problem must 

include at least some v^ -^ Vg oscillations, which will be testable in the near future. 

k>( \ Moreover it also predicts that the lightest super symmetric particle (LSP) decays 

;h ' visibly via lepton number violating modes, which could be searched for at present 

and future accelerators. 



1 Introduction 



It is well-known that within the standard framework with only three light 
neutrinos it is impossible to reconcile solar and atmospheric neutrino data 
[1] with those of LSND [2]. A solution to the atmospheric neutrino problem 
requires a Am^ of the order of Am^ ~ (few)10~^ eV"^ [3,4], while solar neutrino 
data could be explained either by Am^ ~ C(10~^) eV"^ or Am,^ ~ C(10~'^) 
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eV"^ [3,5] (LMA or LOW solutions of the solar neutrino problem), none of 
which can be reconciled with the scale Am"^ ~ 0.3 — 1 eV"^ indicated by the 
LSND experiment [2]. Since the LSND experiment has not been confirmed 
independently, many theoretical papers have chosen to ignore this result 0. 

A simple way to account for all neutrino data would be to postulate the exis- 
tence of an additional neutrino state [7,8], which due to the constraints from 
LEP has to be mainly sterile. 

However, even including a sterile neutrino, there are essentially only two neu- 
trino spectra allowed by the data. One of the possibilities is shown in Fig. (1). 
It consists of two nearly degenerate pairs of neutrinos separated by a gap of the 
order of the LSND scale. The two pairs then have to be split by A?7i? which 
correctly fit solar and atmospheric data. For recent analyses of atmospheric 
and solar neutrino data in four-neutrino models see ref. [9,10]. In Fig. (1) we 
show only the case where the lower Am?- corresponds to the solar scale. A 
second spectrum with Am^^^ <-;> Am^^i is equally well allowed by the data [^ 

The question then is: Can one make sense of such a neutrino spectrum theo- 
retically? 

A number of attempts can be found in the literature [11]. While originally 
motivated by the desire to account for hot dark matter [12-14] soon after 
LSND results came into existence, it was realized that schemes leading to these 
spectra would easily fit LSND results together with solar and atmospheric data. 

Here we propose, what we believe to be one of the simplest particle physics 
model for the inclusion of a sterile neutrino state into the spectrum. It is based 
on a minimal extension of the MSSM with one additional singlet superfield 
field. The fermionic component of the singlet combines with (one of) the ac- 
tive neutrinos to form a light Dirac state at the LSND scale. Its scalar neutrino 
component develops a nonzero vacuum expectation value (vev) [15,16], break- 
ing R-parity spontaneously [17] and effectively generating bilinear R-parity 
violating terms in the superpotential [18,19]. 

Supersymmetry with bilinear R-parity violation has been shown to provide a 
predictive theory for solar and atmospheric neutrino oscillations in which the 
neutrino masses and mixing angles are all determined in terms of the three 
fundamental bilinear terms [20,21]. In the context of the present 4-neutrino 
scheme the breaking of R-parity leads to Majorana masses for the neutrinos, 
splitting the Dirac neutrino into a quasi-Dirac pair [22], and giving mass to one 
additional neutrino state. It is this breaking of R-parity that leads to the solar 
and atmospheric oscillations. The model has therefore, despite being minimal- 



^ The KARMEN experiment [6] does rule out some parts of the parameter space 
favoured by LSND. However, it does not disprove LSND. 

^ The case where Atti^^^ and Am^j^ are exchanged is allowed, but can not be 
realized in our model. 
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Fig. 1. Four-neutrino spectrum fitting solar, atmospheric and LSND data. 

istic, all the basic ingredients for solving solar and atmospheric neutrinos in 
addition to the LSND data. Moreover it extends the predictivity of the bilin- 
ear ^p model to the 4-neutrino case: four parameters are predictable, out of 
a total of ten (4 masses and six mixing angles)[^. In contrast with the early 
models, here the smallness of the LSND scale arises without appealing to a 
radiative mechanism, but due to the volume factor associated with the canon- 
ical normalization of the wave-function of the bulk field in the compactified 
dimensions, as suggested in ref. [23]. 

This paper is organized as follows. In the next section we describe the main 
features of the model. Then, we discuss some numerical results, exploring what 
are the parameter ranges which could fit the neutrino data. We also discuss, 
how this model can be tested by both, future neutrino measurements as well 
as at accelerators, before giving a short conclusion. 



2 The model 



The model we advocate can be regarded as a simple extension of the bilinear 
R-parity broken MSSM [18,19]. The minimal extension beyond the MSSM 
particle content is to include a single right-handed neutrino. Thus to the su- 



We will assume CP invariance throughout this paper. 



perpotential of the MSSM we add, 

W = Wmssm + KU%u\ (1) 

The last term in the equation above gives three Dirac mass terms 

d, = K(vu (2) 



giving automatically a Dirac mass to one linear combination of z/g — z/^ — Vr 
once Hu develops a nonzero vev. The size of this entry is governed by the 
magnitude of the corresponding Yukawa couplings h\. 

However, within a supersymmetric framework the h'^LiHyV'^ term can have 
also other consequences. 

In SUSY models the scalar partner in the v'^ superfield can, depending on the 
parameters of the model, acquire a vacuum expectation value {vr) through 
the usual Higgs mechanism, breaking R-parity spontaneously, as proposed in 
ref. [16]. For this all we need is to have its squared mass negative at the weak 
scale, like in the standard model Higgs mechanism. First note that this may 
happen as a result of the primordial choice of non-universal boundary condi- 
tions at the unification scale which make rn^c different from the SU{2) U{1) 
non-singlet soft scalar masses. Alternatively it might, under some circum- 
stances, be induced dynamically via radiative corrections [24]. 

Such a right-handed sneutrino vev, however, produces effectively bilinear terms 
in the superpotential via 

e^ = Kvr (3) 



It is well-known that the existence of such lepton number violating (R-parity 
breaking) bilinears lead automatically to Majorana masses for the neutrinos. 

In the following we will assume a non-zero vr. With this one can immedi- 
ately write down the mass matrix of the model. Choosing as basis \l'o = 



[iy^..iy„,,iy. 



i\',—i\3,ip}j^,i/j%^) we find that the neutrino-neutralino mass 



matrix of this model is given by: 



Ml ^ ( '- '"«' I . (4) 



where A^yO is the MSSM neutrahno mass matrix, 
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and the Dirac mass matrix tjid is given by, 
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Finally, m^ contains the entries induced by the breaking of R-parity. These 
can be read off from the mass matrix of the spontaneous model as. 
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Here (i/j) are the left-handed sneutrino vevs, which are in general non-zero 
due to the minimization conditions of the scalar potential once we have non- 
zero ej. Note that both va^^ and vhd have such a structure that with either 
term present only one neutrino would gain a mass. Only if mo and m^^ are 
not completely "aligned" to each other (i.e. if at least one of the ratios hijhj 
differs from Aj/Aj where Aj = eiVd + ^{i^i)) we have three non-zero masses plus 
one massless state in the spectrum. 

Depending on the relative size of |Ap/M|^^y and {h^l, where. 



\h' 



lAHy'A^ + A^+A^ 



(9) 



the heaviest state will be either a quasi-Dirac pair or a Majorana neutrino. 
Possible parameter choices which fit the neutrino oscillation data are discussed 
in the next section. 



Before closing this section let us comment on the required smallness of h'^. One 
way to avoid having to simply postulate as a phenomenological assumption is 
to assume the existence of extra compact dimensions, probed only by gravity 
and possibly gauge-singlet fields, which can lower the fundamental scales to 
the weak scale (TeV) [25]. There have been many recent papers applying the 
idea of extra dimensions to neutrino physics [26] . Using the Gauss law one can 
write 

Mil - (R MfTMI, (10) 



where R is the compactification radius of the additional dimensions and Mp 
is the fundamental Planck scale, which in these theories can be low. A recent 
attempt in this direction has been presented in ref. [23]. It postulates that 
Mp — 10 TeV and a 4+n dimensional theory with n = 6 for which the cor- 
responding value of -R is i? ~ 10~^^ cm. This way the smallness of h^ (and 
hence of the LSND mass scale) will follow from the volume factor associated 
with the canonical normalization of the wave-function of the bulk field in the 
compactified dimensions. Our fourth light neutrino Ug {s for sterile) is identi- 
fied with the zero mode of the Kaluza-Klein states. To first approximation the 
sterile neutrino combines with a combination of the active neutrinos (mainly 
i^fj, with some Ur ) in order to form a Dirac neutrino with mass in the eV range 
leaving the other two neutrinos massless. Thus we can apply exactly the same 
construction in the present case. In other words, the present model may be 
regarded as a variant of that in ref. [23] in which neutrino mass splittings are 
now due to the breaking of R-parity. As we discuss below this has important 
phenomenological advantages. 



3 Numerical results 



We now turn to the pattern of neutrino masses and mixings arising from our 
model. Starting with the masses we show in Fig. (2) the eigenvalues which 
follow from diagonalization of the mass matrix for a specific, though arbitrary 
choice of parameters as a function of |A|. For small values of |A| two neutrino 
states form a Dirac pair. Increasing the size of the R-parity breaking (while 
keeping other parameters fixed) increases the mass splitting within this pair, 
as well as the mass of the second mass eigenstate lying at the solar neutrino 
scale. In this example, for |A| ~ 0.05 GeV^ the three mass squared differences 
are of the right order of magnitude for solving the neutrino problems. Note, 
that the actual values of the MSSM parameters are not essential. Larger or 
smaller values of MSSM masses could be accounted for by a simple appropriate 
rescaling of |A|. 

One sees from Fig. (2) that for small values of |A| u^ and u^ combine to form 
a quasi-Dirac pair. The splitting between those two states increases with in- 
creasing |A|. 
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Fig. 2. Example of calculated mass eigenvalues and corresponding mass squared 
differences as a function of |A| for an arbitrary but fixed choice of other model 
parameters. To the left, mass eigenvalues. To the right, Amlg) ^'^Is ^^^ ^"^ii- 
Parameters which have been kept fixed for the calculation shown in this plot are 
Ih^l = 5 X 10-1^ /i = 500 GeV, M2 = 200 GeV and tan/3 = 2.5. 

Turning now to mixings, a general weak interaction four-neutrino model re- 
quires, in addition to the four masses, six mixing angles (and six CP phases 
if CP is violated) in order to characterize the structure of the weak leptonic 
current [8]. 

As explained above the present model has a high degree of predictivity, since 
now all 10 neutrino oscillation parameters are given in terms of six indepen- 
dent quantities which may be conveniently chosen as the di and the alignment 
parameters Aj. Due to the small left-handed sneutrino vevs induced by the 
minimization of the scalar potential the Aj ratios are not proportional to the 
hi ratios, thereby breaking the projectivity of the neutrino mass matrix. This 
is crucial to induce the solar neutrino oscillations. 

Within any model producing a (2, 2)-neutrino spectrum as shown in Fig. (1) 
the conversion probabilities relevant for reactor, atmospheric and LSND exper- 
iments are modified with respect to the usual 2-generation formula. ^ In order 
to compare the experimental data with the results of our model we therefore 
give in the following the conversion (and survival) probabilities corresponding 
to the neutrino spectrum of Fig. (1). 

For the reactor neutrinos one finds after some trivial algebraic manipulations. 



PZi"".?' = 1 -4sin^(AmL^^|)|(l - U^, - Ul){Ul, + Uj,) 



^ Also the solar neutrino survival probability will receive a modification. However, 
CHOOZ and atmospheric neutrino results imply that this correction is negligible. 
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Fig. 3. Two examples of {U^^ + U^^) as a function of /ii//i2 for different choices of 
/i3//i2. This quantity is constrained by CHOOZ [27] and imphes hi < /i2. Only a 
rather weak dependence on other model parameters is found. 
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In Fig. (3) we show two examples of (f/^g + U^^) as a function of /ii//i2 for 
different choices of /i3//i2- This quantity is constrained by the non-observation 
of oscillations at reactor experiments and has to be rather small, implying that 
hi < /i2. 

Similarly for the LSND oscillation probability one finds, 



^^f_^J^ = 4sm\Amls^^-)\{U,iU,i + U,2U,2){U,^U,, + U,,U,,) 
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Since LSND is sensitive only to the largest mass gap in the spectrum, this 
corresponds simply to a re-interpretation of the effective mixing angle. In Fig. 
(4) we show sm^ {29 IS ND,eff) as a function of /ii//i2 for different values of h^lh2. 
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10-^ hi/h2 ~ 0.02 - 0.1 is needed. 



Slightly more complicated is the conversion probability for the atmospheric 
neutrinos, given as 
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Fig. 4. The "effective" LSND angle for the neutrino spectrum of Fig. (1) as a 
function of /11//12 for different choices of /13//12. Obviously to explain the LSND 
result /ii <C /i2 is required. 

The first term proportional to sin^(Am|5^£)^) will be averaged over in the 
atmospheric neutrino data (except for the smallest L and largest E) and 
should be visible in the data as a finite "offset" proportional to 2W{\J^\Ur\ + 

U^2Ur2)iU^3Ur3 + U^4Ur4)\ + \r ^ s|}. The term A{\U^3Ur3U^,4Ur4\ + \r^ s|} 

defines the effective atmospheric neutrino mixing angle in our model. 

In Fig. (5) we give the "effective" atmospheric neutrino angle (left) and the 
corresponding "offset" (right) as a function of h^/h2 for different choices of 
hi/h2. The size of this offset depends mainly on h^/h2- Note that both reactor 
and LSND data require hi <^ h2. 

Note that in the (unrealistically) extreme case of /ii, /is ^ we would have 



pure (two generation) i>^ 



oscillations in our model. However, relatively 



large values of /13 are not excluded by the data, implying that the atmospheric 
neutrino oscillations in our model are described in general by a mixture of 
^tJL — >■ J-'s and Ufj^ ^^ Ut-. Note also that pure u^ — > Ur oscillations are not allowed 
in our model, since the quasi-Dirac pair between which the oscillation takes 
place necessarily involves the sterile state. We have estimated that at least 50 
% of the oscillation probability must be due to z/^ -^ Ug in our model. 

It has been stated that present atmospheric neutrino data rule out a pure 
^'/x -^ ^s conversion [1]. However, we find it premature to confidently rule out 
this possibility at this stage. Indeed, from a global fit in which uncertainties 
are treated conservatively [4], this exclusion does not yet emerge, although the 
pure sterile conversion is indeed disfavored. However, in four-neutrino models 
like the present, the atmospheric conversions are certainly not pure i/^ — *> z/^ in 
general, as seen from eq. (13). As a result, even if ones rules out pure z/^ -^ Ug 
conversion one does not rule out the model itself. However it is definitely clear 
that it will be tested by more refined data yet to come. The bound on the 
admixture of z/„ — > Vg can be inferred from the atmospheric data as suggested 
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Fig. 5. The "effective" atmospheric neutrino angle (left) and "offset" (right) as a 
function of /13//12 for different /11//12 choices. 
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Fig. 6. The "effective" solar angle as a function of Ag/A.^ for two different choices of 

/l3//l2. 

in [9]. Moreover, the "offset" predicted to exist in our model, if the atmospheric 
oscillations are not pure z/^ — !> z/^, might be testable by the K2K experiment. 

To close the discussion on neutrino angles, in Fig. (6) we display the effective 
solar angle as a function of Ag/A,- for two different choices of other parameters. 
Due to the flatness of the recoil electron neutrino spectrum indicated by the 
most recent solar neutrino data [1], the solutions to the solar neutrino problem 
which are now preferred by the data involve large mixing [3,5]. Large neutrino 
mixing angles require Ag ~ A,- within a factor of 2 — 3. 
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4 Neutralino decay 



Let us finally briefly discuss the decay of the lightest of the 4 heavy states 
in our neutrino-neutralino spectrum. We will call it for simplicity the lightest 
neutralino. As is well-known with the presence of the bilinear (lepton number 
violating) terms in the superpotential the lightest neutralino is no longer stable 
[28,29] and decays through lepton number violating modes. In the bilinear R- 
parity violating model it has been shown that this decay occurs with sizeable 
branching ratio into visible states [20] inside typical detector sizes despite 
the small neutrino masses required by current neutrino data. Moreover, these 
decays might be used to obtain information on neutrino angles [30]. 

Since in our present model R-parity violation is required in order to provide 
the solar and atmospheric neutrino mass splittings, the neutralino will decay 
as well. However, two features distinguish the current model from the simpler 
bilinear model. First, the spontaneous violation of R-parity (which is used as a 
seed for the generation of the bilinear terms) implies the existence of a Majoron. 
The neutralino then also can decay via the invisible mode x° -^ Jv [28], see 
Fig. (7). 

However, due to the smallness of the (Dirac) neutrino coupling this decay 
mode will be very much suppressed when compared to the visible ^p violating 
decays. A very crude order-of-magnitude estimate for the two graphs shown in 
Fig. (7) gives for the neutral current graph 

r- ^ ^(^)2^ ,, ^ 10-16 GeV (14) 

while for the Majoron decay one expects: 



m^o ~ 10"'' GeV (15) 



IGvr 

Thus, even though a Majoron exists in the present model, sufficiently large 
branching ratios to visible states in the decay of x" exist to be searched for at 
accelerators. 

Second and more important, however, in the present model the solar neutrino 
problem is solved by z/g -^ v^ oscillations (with possibly some admixture of 
^e ~^ '^n)- If the solar neutrinos indeed require large mixing, as currently 
favoured by the data, Ae — A,- is needed in our model to fit the data. This 
would then lead to the prediction that in the decays of the neutralino to 
semi-leptonic final states a comparable number of electrons and taus should 
exist. This prediction might be used to distinguish the present model from the 
bilinear R-parity violating model, since the latter requires that electrons in the 
semi-leptonic states should be very much suppressed [20,30] refiecting the fact 
that in the bilinear model Ae ^ A^ ^ A,- is required. 
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Fig. 7. Two examples of Feynman graphs leading to the decay of the lightest neu- 
tralino in our model. To the left, neutral current interaction leading to visible final 
states. To the right, invisible decay into a Majoron and a neutrino. 

5 Conclusion 



We have discussed a very economical model of neutrino mass which, despite 
its simplicity is in principle able to fit all neutrino data including LSND. 
The model is minimal in the sense that we only introduce one additional 
SU{2) (g) f/(l) singlet superfield into the MSSM superpotential. The LSND 
scale can be explained by the fermion present in this superfield, while its scalar 
partner can trigger the spontaneous violation of R-parity, through a nonzero 
right-handed sneutrino vacuum expectation value, which in turn splits the 
Dirac neutrino into a quasi-Dirac pair (atmospheric neutrino scale), provid- 
ing also mass to one of the low-lying neutrinos at the solar neutrino scale. 
The smallness of the overall neutrino mass (LSND scale) can be ascribed to 
the volume suppression factor associated with the compactified dimensions. 
Thus the model can naturally explain all neutrino oscillation data, making 
four predictions for the neutrino oscillation parameters. We stress that this 
model is testable in the near future: one of its necessary ingredients is that 
atmospheric neutrino conversions must include at least some v^ —>■ Vg oscilla- 
tions. Thus future atmospheric neutrino measurements will provide a crucial 
test of this model. However, if tides turn, and sterile neutrinos in the future 
are not disfavored by Superkamiokande, another test of the model could be 
done at accelerators. The model predicts that the lightest neutralino decays 
inside the detector with sizeable branching ratios into visible states. 
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